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t— I ■ Abstract 

^ ; In the SND experiment at VEPP-2M e+e" collider the 0(1020) -> 7r°7r°7 decay was 

. studied and its branching ratio was measured: B {4> -» 7r°7r°7) = (1.221±0.098±0.061)-10- 4 . 

<0 It was shown, that f (980)7 intermediate state dominates in this decay and the /o(980)- 

meson parameters were obtained. 
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j>! ; 1 Introduction 

•rH , 

The first work to search for the decay 4> — > 7r 7r°7 was performed in 1987 |l], ^|. The observation 
of this decay based on analysis of the data of the PHI96 experiment Q was reported in 1997 
|Q. |5[ by the SND detector §. Observed decay rate of about 10 -4 was unexpectedly high and 
it was also shown that /o (980)7 intermediate state dominates in this decay. Later these results 
were confirmed by the CMD-2 experiment . 

In this work results are presented of a new experimental study of the process 

e+e" -» ^ -» 7r°7r°7 (1) 

based on complete SND statistics of about 2 • 10 7 produced 4> mesons. The experiments were 
carried out at VEPP-2M e + e~ collider in Novosibirsk. The total integrated luminosity collected 
in PHI96 ||] and PHI98 Q experiments was 4 pb' 1 in 20 energy points and 8 pb 1 in the 16 
energy points respectively. 

*The work is partially supported by RFBR (Grants No 99-02-16813, 99-02-16815, 00-02-17478, 00-02-17481) 
and "Russian Universities" (Grant No 3H-339-00). 
' email : V . N . I vanchenko® inp .nsk.su 
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2 Data analysis 



Main resonant background to the reaction (|l|) comes from the process 

e + e~ — ► <f) — ► 777 — ► 37r°7 (2) 

due to the merging of photons and/or loss of photons through the openings in the calorimeter. 
A smaller background comes from the decay mode || 

e + e~ — > — > r/7r°7 (3) 

Copiously produced </> — ► K$Kl —* ^ir^K^ events may contribute to the background but this 
decay is almost completely rejected by our cuts. The main source of non-resonant background 
is the process 



e e — ► W7r 



. „o_o 



7r U 7r u 7 . (4) 

The QED background from e + e~ annihilation into 4,57 was found negligible. 

The data analysis was described in detail in our previous publications 0, |5|. In this work 
basically the same selection cuts were used requiring five isolated photons in the calorimeter 
with good shower profiles and event kinematics consistent with the hypothesis (93 % 

c.L). To suppress spurious photons from beams background all photons with energies less than 
60 MeV at angles with respect to the beams less than 36 degrees were excluded from further 
analysis. The rate of events with spurious photons was about 4 % in PHI96 experiment and 8 % 
in PHI98. To account for this effect during processing of simulated events spurious hits, taken 
from the experiment, were admixed to Monte Carlo ones. 

After applying cuts the sample of the 7r°7r°7 candidates was selected (Fig.|l]a). The mam 
background is due to the process (0) but processes (Q) and (|3|) contribute to this event sample 
also. The estimated number of events of the reaction ([l]) in FigJTja is 419 ± 31. The detection 
efficiency (Fig.[l]b) for the process (|l|) as a function of 7r°7r° invariant mass was determined using 
simulation of the process <fi — > where 5 is a scalar state with a mass ranging from 

300 to 1000 MeV. The whole mass range was divided into 19 bins and within each bin simulation 
with uniform distribution was performed. This simulation was used for determination of 
the 7r°7r° invariant mass resolution and event misidentification probabilities as well. 

For the background process (Q) the 7r°7 invariant mass distribution is peaked at w-meson 
mass. The m 7r7 parameter is an invariant mass of the recoil photon and one of 7r° mesons, closest 
to the w-meson mass (Fig.||). Events with 

750 MeV < nvy < 815 MeV (5) 

are mainly ones of the process (H) . Moreover in the energy range of the experiment vr 7r° invariant 
mass in the process (||) must be lower than 700 MeV. For further analysis the selected 7r°7r°7 
events were divided into three classes: 

1. < 700 MeV and not (|); 

2. m n7T > 700 MeV and not dBD; 
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3. m W7T < 700 MeV and (§). 

In the second class the /07 mechanism dominates, in the third — the process (Q). 

For events of the second and third classes the experimental and simulated distributions in 
tp and angles were compared. The ip is an angle of the recoil photon with respect to pion 
direction in the 7r°7r° center of mass reference frame. The 9 is an angle between recoil photon 
and the beam. It can be seen in Fig.||a,c that in the second class the distribution well matches 
one expected for the intermediate scalar state. Corresponding distributions in the third class 
agree with ones expected for the intermediate ujtt state (Fig]3|b,d). 



3 Fitting of cross sections 

The visible cross sections (Ji(s) were measured for each energy point in both PHI96 and PHI98 
experiments, here s = 4E 2 and i is a class number (i=l,2,3) defined above. These cross sections 
were fitted according to the following formulae: 

00, (6) 

3=1,2 

where er£ (s) is a cross section for the j-th m wn interval (J = 1 for m n7T < 700 MeV, j = 2 
for > 700 MeV), are probabilities for the 7r°7r°7 events belonging to j-th m w ^ interval 
to be detected as the class i events (Tabl|l|,§), 07 (s) are production cross sections of the final 
state / (/ = 777, ?7vr°7, lott ), are detection efficiencies of the state / in the class i, a/(s) - 
radiative corrections and corrections to the energy spread in the beam |10|. The expression for 



the (Jrn{s) was taken from Ref.[ll], the values of the a LU7T (s) as well as a UJ7T (s) were taken from our 
previous work [12]. The efficiencies e l ^ n ^ were estimated using the experimental mass spectrum 
and the simulation of £7 events described in the previous section. For the cross sections 07 (s) 
(f is 7r°7r°7 or 7r°r]j) the following expression was used: 

127^5(0 - ee)B{<j> -> Rj)m%F(8) 
a f (s) = v — 



F(s) = ^ R fL R = fo,ao- 



Here the inverse propagator of was written according to Ref.[13]: D^s) = m| — s — i- y / ir^ ) (s) 
and the energy dependence of partial widths T^r^^s) was taken from Ref.[14]. The branching 
ratio of the decay (||) was obtained in SND work || . 

Note, that there is a dip in the cross section (|) at the (i-meson mass, due to destructive 
interference between the main non-resonance process, the </> — ► lutt° decay, and the process 



e + e~ -> <f> -> /9vr u -> vr u 7r u 7. (7) 

To check the accuracy of the detection efficiency, obtained by simulation, fully reconstructed 
events of the process @ with 7 photons in the final state were analyzed. Since this process is 
similar to (ffl), the energy, invariant mass, and angular cuts were the same as in 7r 7r°7 analysis, 
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while threshold values for kinematic fitting parameters were scaled to provide same confidence 
level. Obtained branching ratio of the decay (|2|) was smaller than that from the SND experiment 



B(<f> 


-» ttVt) = 


(1.034 ± 0.066 ± 0.046) 


io- 4 , 


m n7T > 700 MeV, 


B(<f> 


-» vr 7r°7) = 


(0.559 ±0.053 ±0.025) 


10~ 4 , 


> 900 MeV, 


B(4> 


- 7r 7r° 7 ) = 


(0.124 ±0.065 ±0.006) 


• io~ 4 , 


< 700 MeV. 



lTJ B(4> -> 7/7) = (1.338 ± 0.053)% by a factor of £ = 0.89 ± 0.04. This value was taken as 
a correction factor for both effect and background. This correction accounts for imprecise 
simulation of tails of distributions of parameters used in the cuts. 

Free parameters of the fit were the branching ratios in the 0-meson peak B((f) — > 7r°7r°7) 
in two intervals of m W7V (Tabljl],^). Uncertainties in the values of r^, detection efficiencies, 
integrated luminosity, and energy spread of the beam were taken into account. The value of 
B((f) — ► 7/7) was fixed to its SND value, so our results are normalized to the branching ratio of the 
decay <fi — ► 7/7. The fitting results for both experiments are shown in the Table ||, corresponding 
curves for PHI98 are shown in Fig.|I[ Because the most of 7r°7r°7 events concentrate at high 
similar analysis was performed for the mass range > 900 MeV. 

4 Results 

Because the fitting results for the PHI96 and PHI98 experiments are in a good agreement with 
each other (Tabl. |||) they was combined: 

(8) 
(9) 
(10) 

In (|8]-4To|) the first error is statistical and the second one is systematic, which is determined 
mainly by current accuracy of B((f> — > 7/7) (4 %) and model uncertainty (2 %) of F(s) in the 
formula (g). From the data @ and @ one can evaluate the total value 

B(<j) -> 7r vr°7) = (1.158 ± 0.093 ± 0.052) • 10~ 4 , -> wvr excluded, (11) 

which is the total branching ratio excluding — ► a;7r° contribution, because the kinematic region 
of this decay is excluded by our cuts. It was studied in the SND work from which one can 
obtain the branching ratio to the unr region B(<j) 7r°7r°7) = (0.063 ± 0.031) • 10 -4 and add it 
to dTI| ) to find the total branching ratio 

B((f> -» vr vr 7) = (1.221 ± 0.098 ± 0.061) • 10" 4 . (12) 

In this work the branching ratios were measured relatively to B(<p — * 7/7). In such approach 
some systematic errors cancel allowing smaller total systematic uncertainties: 

(0.913 ± 0.073 ± 0.029) • 10~ 2 , total; 
(0.865 ± 0.070 ± 0.017) • 10~ 2 , ujtt excluded; 
(0.773 ± 0.049 ± 0.016) • 10~ 2 , > 700 MeV; 
(0.418 ± 0.040 ± 0.009) • 10~ 2 , > 900 MeV. (13) 



B(cf>- 


> 7T 7r°7) 


B(cf> 


-» 7/7) 


B{4>- 


> 7T 7T 7) 


B{<t> 


-» 7/7) 


B{4>- 


> ^°7T 7) 


B{<t> 


-» 7/7) 


B(cf>- 


> ^°7T 7) 


B(cf> 


-» 7/7) 
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To investigate the dynamics of the reaction (|T]) the 7r°7r° invariant mass spectrum (Fig.[l|a) 
was studied in a narrower energy interval 2E = (1015 — 1025) MeV containing almost all cp- 
meson events and smaller ujtt background. The estimated background from the processes (|| [|) 
was subtracted using SND data g H, Q. Resulting raw mass spectrum is distorted due to 



finite detector resolution and misidentification of recoil photon in some 7r°7r°7 events. In order 
to obtain corrected spectrum the matrix ej j was constructed, where each ej j is a probability to 
get a reconstructed value within the invariant mass bin i for an simulated event with actual 
m n7T uniformly distributed in the bin j. Note, that in the two highest bins the simulated mass 
distribution were subdivided into smaller bins to take into account rapid decrease of the mass 
spectrum. The diagonal elements of the matrix and their immediate neighbors describe detector 
resolution. They turned out to be much larger than others which are due to misidentification of 
the recoil photon. Using this matrix and the experimental mass spectrum the number of events 
with misidentification were estimated and subtracted from each bin. This procedure does not 
affect the diagonal elements of the matrix and their immediate neighbors but zeroes all others. 

After normalization to the total integrated luminosity and detection efficiency the differential 
branching ratios S exp (mi) = dB{4> —* "f)/dm n7 r for both experiments were obtained. They 
are in a good agreement and thus can be combined together (Tabl.||). Using the matrix of 
efficiencies it is possible to restore the physical mass spectrum but such a procedure makes 
the values of corrected spectrum strongly correlated. To compare these results with model 
predictions one has to deal with a cumbersome covariance matrix. To simplify the analysis we 
instead transformed bare theory predictions Sth{ m ) using the values of efficiencies from Table |j: 

Scorim) = €i,i-iS t h(mi-i) + e it iSth(mi) + ei, i+1 S th (rn i+ i), (14) 

here S cor (rrii) is a corrected for the detector response theoretical mass spectrum, which should 
be used for approximation of the experimental spectrum S exp (mi) by the least squares method. 
This approach is used in our fits, which are discussed in the next section. 

5 Discussion 

The 7r°7r° invariant mass spectrum (Figj^) was approximated by the following function: 



Sth(m) = Ssryim) + S pn (m) + 2^J S Sl {m)S p7 ,{m) cos 0, (15) 

where <S , 5 7 (m) is the contribution of the /o7 and o"7 decays, S pn (m) is the contribution of the 
process ([?]), and is an interference phase. This formula is exact only in case of identical 
angular distributions in Sj and 

A 

intermediate states. Provided small contribution, the 
difference in the angular distributions for the 7r°7r°7 final state can be approximately accounted 
for by fitting of 9. The S pn (m) dependence was taken from [14|, [T(|. Corresponding branching 
ratio is B(<j) pn vr 7r°7) = (1.2 ± 0.2) • 10 -5 |]|. For S s ~/(m) the expression from Ref.||] 
was used 



Ss~[(m) 



2m 2 T 



4>f7 L /7r°7r° 



1 A 

+ 



Df(m) D a (m) 
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where Ds(rn) = m 2 s — m? + lis (to) — i ■ mTs{m) is /q or a inverse propagator, lis (to) takes into 
account finite width corrections |fL4]| , A is a complex parameter taking into account difference in 
coupling constants of /o and a mesons. 



In the kaon loop model of the <fi — ► /07 transition [jlj] r<^/ 7 is proportional to the product 
g'j K+K _g'j n+7T _. The dependence of r^/ 7 on to in this model is different from the model with 
the point-like — ► /07 transition and can be interpreted as a form factor, which stops the r^/^ 
growth at higher photon energies. 

The free parameters of the fit were 6, A, nif, and coupling constants gjx+K- an d gfn+n- 
while the width and mass of the a meson were fixed at m a = 600 MeV and T a = 400 MeV (TtJ. 
For a meson the Breit-Wigner propagator with II -(to) = and T a (m) = const was used. The 
fitting result is A = — 0.065±g;?|^, indicating that aj contribution is compatible with zero. This 
allows to reduce the number of free fitting parameters by fixing A = 0. For the remaining free 
parameters the following values were obtained: 

to/ = (969.8 ± 4.5) MeV, 
g} K+K -/A7r = 2.471°;™ GeV 2 , 
g% +n -/4ir = 0.54i^GeV 2 , 

9}k + k-/9% +7T - = 4.6 ±0.8, 

9 = (180 ± 36) degrees, 
X 2 /n D = 3/14. (17) 

Note, that although contribution from the process (0) is small, and fitting with this contribution 
set to zero still gives high confidence level (x 2 / n D = 10/14), but corresponding fitting curve 
(fig.||) shows visible systematic deviations from the observed mass spectrum. 

The fitting results (|l7|) demonstrate that the _fo7 mechanism dominates in the decay (Q) , and 
contributions from (17 and p°Tr° are small. Neglecting these contributions and assuming natural 
isotopic ratio B(fo — > -K + ir~) = 2B(fo — ► 7r°7r°) one can obtain from ( JTT| ) 

5(0^/o7) = (3.5±0.3i^)-10- 4 , 

W~r|| = (2-6±0.2l°|)-10- 2 . (18) 

where the second errors are systematic which includes uncertainty in the interference terms. 
These errors were estimated by summation of systematic uncertainty from (O) and the difference 



between (18) and the value B((f> — > /07) = 4.6 • 10 4 calculated from (1171 ) with (17 and 
contributions taken into account. 

The /o- meson width in the model ||T3 is a rapid function of to. At the resonance pole 
it is equal to Tf(mf) = —Im(Df(mf))/mf = (201 db 28) MeV. Due to its close proximity 
to KK threshold this value is significantly larger than visible in Fig.|| /o width ~ 100 MeV. 
Note, that if the width corrections to the real part of /q propagator are set to zero the values 
of fitting parameters including /o mass and width change significantly: to/ = 985^^ MeV, 
Tf = 122 ± 13 MeV, but such assumption looks much less attractive from the theoretical point 
of view. The model dependence of the /o mass and width explains disagreement between (fij 



and some data from PDG [18]. 
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The model of Ref . |14| which takes into account the kaon loop contribution in the decay under 
study fits well the experimental spectrum. On the contrary, it is not possible to describe the 
data using the model with point-like transition mechanism, in which the transition form factor 
is equal to unity and T(j,f -y is proportional to w 3 , where u> is the photon energy. In this case 
the fitting curve does not follow the experimental spectrum and the x 2 / n D = 28/14 (fig.|5|). So, 
kaon loop contribution, which provides non point-like transition (j) — > /07 is confirmed by the 
our data. 

The values of coupling constants ( |17| ) are in a good agreement with the 4-quark MIT bag 
model predictions 0, ||, |(|. New theoretical studies |l|, H, ||, ^ |j| were performed after 



first observation in order to explain the data. It was shown that it is possible to describe exper- 
imental results using different theoretical approaches. In these models the 4-quark component 
exists directly or "de facto" as a result of the kaon loops contribution. 

6 Conclusion 

Results of this work confirming our first publications @, |f| and CMD-2 measurements [Q] are 
based on the full SND data set which contains the largest statistics of 7r°7r 7 events. The 
shape of the 7r°7r° mass spectrum is described well by the sum of pir and /07 contributions. 
The contribution of a'j intermediate state is not required for description of the data but is not 
completely excluded by our data. It is shown that point-like transition model does not fit the 
experimental mass spectrum. At the same time, models with intermediate kaon loops reproduce 
the data. Predictions of the four-quark MIT bag model of the /o meson are in a good agreement 
with our results but we cannot exclude other explanations. Further theoretical studies and new 
data are needed to obtain consistent description of the /o(980)-meson quark structure and to 
detect o"7 contribution in the (f) — > decay. 
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Table 1: The matrix of detection efficiencies for three classes of events for the PHI96 experiment. 



3 


1 


2 


m^, MeV 


300-700 


700-1000 


e lj % 

c 7T7r7> /u 

e 2j % 

^7T7r7 1 ' u 

e 3j % 

7T7T7 ' U 


13.27 

1.46 

5.11 


0.82 

18.25 

0.45 



Table 2: The matrix of detection efficiencies for three classes of events for the PHI98 experiment. 



3 


1 


2 


m nn , MeV 


300-700 


700-1000 


C 7T7T7> /U 

e 2 i % 

u 7T7r7' /u 

e 3j % 

7T7T7 ' U 


14.34 


0.87 


1.79 


18.87 


4.19 


0.47 



Table 3: The results of the fit of the cross sections. 



Experiment 








PHI96 


PHI98 


B(</> 7T U 7T U 7) 


•10 4 , 




< 700MeV 


0.258 ±0.111 


0.055 ± 0.080 


-► vr 7r°7) ■ 


• 10 4 , 


fn n7T 


> 700MeV 


0.960 ±0.102 


1.087 ±0.086 


B{4> -► 7r°7r°7) 


■ 10 4 , 


fn n7T 


> 900MeV 


0.582 ±0.085 


0.545 ± 0.068 
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Table 4: The m wn mass bins, average mass inside the bin nii, the non zero elements of the 
efficiency matrix eij, and the normalized experimental invariant mass spectrum for the <fi — > 
7r°7r°7 decay after background subtraction and acceptance corrections. Only statistical error are 
shown. Systematic uncertainty is 4.5 %. 



N 


m n7T , MeV 


rrii , MeV 






e i+l,i 


dB V m ir i] ■ 10 8 (MeV- 1 ) 


1 


300-400 


350 




0.953 


0.047 


3.36 ± 4.47 


2 


400-500 


450 


0.037 


0.920 


0.043 


1.01 ±4.67 


3 


500-600 


550 


0.043 


0.900 


0.056 


1.94 ±3.11 


4 


600-700 


650 


0.042 


0.718 


0.240 


6.26 ± 2.00 


5 


700-720 


710 


0.058 


0.660 


0.281 


9.46 ± 6.46 


6 


720-740 


730 


0.227 


0.536 


0.237 


13.99 ±6.14 


7 


740-760 


750 


0.221 


0.544 


0.235 


13.53 ±6.07 


8 


760-780 


770 


0.240 


0.537 


0.223 


18.23 ±6.26 


9 


780-800 


790 


0.242 


0.546 


0.212 


16.93 ± 6.08 


10 


800-820 


810 


0.221 


0.558 


0.220 


24.04 ± 7.36 


11 


820-840 


830 


0.233 


0.559 


0.208 


28.07 ±7.65 


12 


840-860 


850 


0.223 


0.572 


0.205 


26.59 ± 7.59 


13 


860-880 


870 


0.215 


0.564 


0.221 


34.15 ±8.28 


14 


880-900 


890 


0.202 


0.613 


0.185 


44.37 ± 9.37 


15 


900-920 


910 


0.192 


0.630 


0.179 


43.15 ±8.95 


16 


920-940 


930 


0.197 


0.651 


0.153 


61.66 ± 10.38 


17 


940-960 


950 


0.171 


0.698 


0.131 


65.45 ± 10.71 


18 


960-980 


969 


0.178 


0.728 


0.095 


70.14 ±12.61 


19 


980-1000 


987 


0.214 


0.786 




30.52 ± 10.08 
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Figure 1: The mass spectra for events: a — the invariant mass distribution of 7r°7r° 

pairs without acceptance corrections, beam energy 1015 MeV < 2E < 1025 MeV, histogram - 
data, shaded histogram — estimated background; b — the detection efficiency for the process 
e + e~ — ► Sj — > 7r°7r°7, S is a scalar state. 
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900 



m nr MeV 



Figure 2: Distribution of 7r°7 invariant mass for < 700 MeV. Points — data, histogram 



simulation, shaded histogram 



sum of contributions from the <f> decays, arrows — selection of 
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0, degrees 0, degrees 



Figure 3: The angular spectra for ?r 71" 7 events: a, b — distribution of cosine of ip, the angle 
between directions of tt° and recoil 7 in the rest frame of system; c, d — distributions of 9, 
angle of the recoil 7 with respect to the beam. Points — data, histogram — simulation, shaded 
histogram — sum of contributions from background processes for the reaction (|l|) (a,c) and for 
(U) (b,d). Acceptance corrections are not applied. Beam energy is 1015 MeV < 2E < 1025 MeV. 
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Figure 4: Energy dependence of the visible e + e~ — ► 7r°7r°7 cross section for the three classes of 
events described in the text for the PHI98 experiment. Points — data, solid line — fit, dotted 
line — sum of background contributions to the decay <f> — ► 7r 7r° 7 . 
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Figure 5: The 7r tt invariant mass spectrum after background subtraction and acceptance 
corrections: points — data, solid line - the result of the fit, dashed line — the result of the fit 
without contribution; dotted line — the result of the fit in the point-like transition model. 
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